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Abstract 

A signal of two leptons and missing energy is challenging to analyze at the Large 
Hadron Collider (LHC) since it offers only few kinematical handles. This signature 
generally arises from pair production of heavy charged particles which each decay 
into a lepton and a weakly interacting stable particle. Here this class of processes is 
analyzed with minimal model assumptions by considering all possible combinations of 
spin 0, g or 1, and of weak iso-singlets, -doublets or -triplets for the new particles. 
Adding to existing work on mass and spin measurements, two new variables for spin 
determination and an asymmetry for the determination of the couplings of the new 
particles are introduced. It is shown that these observables allow one to independently 
determine the spin and the couplings of the new particles, except for a few cases that 
turn out to be indistinguishable at the LHC. These findings are corroborated by results 
of an alternative analysis strategy based on an automated likelihood test. 



1 Introduction 



Many models beyond the standard model (SM) include stable weekly interacting massive 
particles which could be constituents of dark matter. Since the stability of these particles is 
generally related to some symmetry, they can be produced only in pairs at colliders, leading 
to challenging signatures with at least two invisible objects. At hadron accelerators like the 
Large Hadron Collider (LHC) such a signal is not sufficiently kinematically constrained to 
enable the use of direct reconstruction techniques, and thus it is very difficult to uniquely 
determine the properties of the produced particles. 

One of the most challenging cases are processes with a low-multiplicity final state of 
only two visible objects, which is focus of this article. In particular we will consider the 
production of a pair of oppositely charged heavy new particles Y ± at the LHC, which each 
decay into a SM lepton and an invisible neutral massive particle X°: 

pp -> Y + Y- ->■ £ + rX°X°, (£ = e,fi). (1) 

Several methods for determining the Y and X masses in processes of this type have been 
proposed in the literature [IHI]- Furthermore, a number of authors have studied how to 
extract spin information from angular distributions [SHI] and the total production rate |8jE 
To the best knowledge of the authors, however, the problem of determining the couplings of 
the new particles, which are related to their gauge group representations, has not yet been 
considered. 

The goal of this article is to analyze the process ([I]) in a more model-independent approach 
by considering all possible assignments for spins (up to spin one) and SU(2) representations 
(up to triplets) for the particles X and Y. To discriminate between these template model 
combinations, we discuss several variables for the measurement of masses, spins, and in- 
teractions of the new particles, including two new spin-sensitive observables and one new 
observable for the coupling determination. To minimize model dependence, the total cross 
section is not considered in this set of variables. 

Besides using dedicated observables, we also study an alternative analysis strategy based 
on an automated likelihood test. This method matches the observed lepton momenta in a 
sample of events to the corresponding momenta of a theoretically calculated matrix element 
in a given model and calculates a likelihood from that 0,011] ■ Note that the two approaches 
based on specific observables and on the automated likelihood test are complementary. The 
latter often reaches a higher statistical significance due to the fact that no information is 
lost by projecting onto some variables, but it does not permit a straightforward separation 
between individual particle properties, such as spin and couplings. 

After characterizing all relevant spin and coupling representations in section [2] and iden- 
tifying 11 independent model combinations, we present observables for the measurement of 
particle properties in section [3] and demonstrate their usefulness in a Monte-Carlo study. 
Section H] is devoted to the analysis of the same set of template models with the likelihood 
test method. Finally, our conclusions are given in section [5j 

*The latter method, while potentially very powerful, requires knowledge of the branching fractions of Y^, 
which are a priori unknown without model assumptions. 
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Figure 1: Basic diagram topology for the new physics processes under consideration. Thick 
lines indicate new particles, while thin lines denote SM particles. 



2 Setup 

The class of processes under consideration each involve Drell-Yan-type production of a pair 
of charged heavy particles Y ± , which each subsequently decay into a lepton £ ± (£ = e, fi) and 
a neutral heavy particle X°, see Fig. [TJ It is assumed that the and X° are charged under 
some discrete symmetry, such that they can be produced only in pairs and the lighter new 
particle (X°) is stable and escapes from the detector without leaving a signal. The observable 
signature thus consists of two same-flavor opposite-sign leptons and missing momentum: 
£ + £~ + ft. For this process it is insubstantial whether X° is self-conjugate or not. 

For the purpose of this work it is assumed that no other new heavy particles play a role 
in the s- or t-channel of the Y + Y~ production process [note, however, that if Y is a vector 
boson this assumption is not valid, as will be explained below]. In fact, LHC data itself will 
be able to set a strong lower bound on such particles: searches for di-jet resonances could 
rule out s-channel resonances that couple to light quarks up to several TeV, and any new 
particles in the t-channel need to be colored and thus could be produced directly with large 
cross sections unless their masses are larger than 2-3 TeV [llj^l . Therefore, if the LHC does 
not see any such signals, one can safely neglect the presence of extra particles in the s- and 
t-channel for the production of Y + Y~ pairs with mass of a few hundred GeV. 

Table [1] lists 16 possible combinations of spins up to spin one and singlet, doublet, or 
(adjoint) triplet representations under the weak SU(2) for the fields X and Y . We do not 
consider complex SU(2) triplets, since they contain doubly charged particles, which would 
lead to a clearly distinguishable signature. 

Also shown are the structure of the couplings between the Y and the Z boson and 
between X, Y, and SM charged leptons. The jYY couplings has the same form as the 
ZYY coupling. The coupling constants for the ZYY coupling are shown in Table [2, given 
in terms of the ratio with respect to the / -/YY coupling, Rza = g(ZYY)/e. The strength 
of the XYi coupling depends on the detail of the given model, and it is only relevant for 
the overall decay branchings, but not for the shapes of distributions. We neglect corrections 
from electroweak symmetry breaking to the masses and interactions of X and Y . As a result, 
if Y is a spin-1/2 fermion it couples to the Z boson only through non-chiral vector couplings. 

For illustration, Tab. [T] also gives examples for concrete realizations of all 16 spin and 



t These estimated bounds pertain to the LHC with a center-of-mass energy of 14 TeV. 
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Table 1: List of different assignments of spin s and SU(2) representations for the charged 
held Y and the neutral held X. We dehne Y~ /Y + to he the particle/ anti-particle. Also 
shown are the structure of the couplings to the Z boson and to SM leptons, as well as 
examples for realizations of these assignments in known models. MSSM refers to the Minimal 
Supersymmetric Standard Model, UED to (at least) one universal extra dimension, and 
UED to (at least) two universal extra dimensions, t^, i^, v, B°, W 0)± , and H denote the 
superpartners of the right-handed charged lepton, left-handed charged lepton, neutrino, U(l ) 
gauge held, SU(2) gauge helds, and Higgs boson, respectively, ig^y i~ D , x y V{i), B° W^'/t, 

and H(i), respectively, are the hrst-level KK-excitations of these helds. B° H ,^ and ^ are 
scalars stemming from one of the extra components of the higher- dimensional gauge helds 
in UED. More details of these models can be found in Refs. fI2 |, [T3 y . 
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Table 2: Ratio Rza of the ZYY to ^YY coupling strength for different SU(2) representa- 
tions ofY. 

gauge group assignments within the Minimal Supersymmetric Standard Model (MSSM) or 
models with universal extra dimensions. However, many of these combinations can also be 
realized in other models. 

A comment is on order regarding the combination 11 in the table. Taking only the s- 
channel diagrams in Fig.[H the cross section for spin-1 Y pair production grows unboundedly 
for increasing partonic center-of-mass energy This is a result of incomplete SU(2) gauge 
cancellations. Gauge invariance requires the presence of an additional new particle in the 
t-channel, which interferes negatively with the s-channel contribution and thus preserves 
perturbative unitarity. In the case of universal extra dimensions this role is played by the 
KK-quarks. Therefore, for model 11 we include a new colored fermion Q that is charged 
under the same discrete symmetry as X and Y. The coupling strength of the qQY interaction 
is prescribed by gauge invariance: g(qQY) = g. While for consistency it is necessary to 
incorporate this particle in the cross-section calculation, it is still possible that it is too 
massive to be seen directly at the LHC, i. e. tjIq > 0(TeV). 

The cross sections for the Drell-Yan-type process in Fig. [1] for the different models in 
Tab. [1] range from a few fb to several hundred fb for a center-of-mass energy of \fs = 14 TeV 
and X/Y masses of a few hundred GeV, see appendix. Thus one can expect several 100- 
10,000 events being produced with a total luminosity of 100 fb -1 . Note that in this paper, as 
mentioned in the introduction, the total event rate will not be used to discriminate between 
models, since it would require knowledge of the Y ± branching fractions. 

At the LHC it is not possible to determine the polarizations of the final-state leptons 
and X particles. As a result, several pairs of combinations in Tab. [1] are indistinguishable 
from each other, since after summing over the spins of the external legs their squared matrix 
elements are identical. Those sets of look-alikes are (1, la), (2, 2a, 2b), (7, 7a), and (9, 9a). 
This leaves a total of 11 potentially distinguishable combinations, which will be explored in 
more detail in the following. 

These 11 combinations have been implemented into CompHEP 4.5.1 [Hj and repre- 
sentative samples with a few thousand parton-level events have been generated for each of 
them. Since we will not consider the total cross section as a discriminative quantity, the 
exact values for the XY£ coupling strength and the widths of the Y particles are irrelevant. 
However, the Y widths have been chosen small enough so that diagrams with off-shell Y 
particles can be safely neglected, i. e. T Y /M Y <C 1%. 

As a first step, initial-state radiation and detector acceptance effects have been ignored 
in the following analysis, but we discuss these contributions in section 13.51 
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3 Observables for Determination of Particle Properties 



The experimental information in the signature £ + £~ + E consists of the 3-momenta of the 
leptons £ + and t~ , which can be parametrized in terms of their transverse momentum pr, 
pseudorapidity 77 and azimuthal angle 0. Since the system is invariant under overall az- 
imuthal rotations, one can construct five independent non-trivial observables from this data. 
In this work we will focus on the following five quantities: 

Mt2 = min jmax (m^ ,Xl ,m e T ' X2 ) i , (2) 
Pt,x 1 +Pt,x 2 =F*t L J 

cos 9* u = tanh ^ + - ^ £ ■ = tanh — ^— , (3) 

M eS = p T/ + + p T/ - + // T , (4) 

A(f) U = \(<f>£+ - <j>f) mod 2tt|, (5) 
. N(E e - > E e +) — N(Ej+ > E e -) 

e+e ~ N(E e - > E e +) + N(E e + > Eg-)' U 



Here 



(rn T ' Xi ) 2 = m 2 x + 2(p T/ ±Jm 2 x + p| Xi - p T/ ± ■ p TX .) 



is the transverse mass of the lepton assumed to be massless, and one neutral heavy 
particle X^, i = 1,2. Furthermore, 9* u is the polar angle between one lepton and the beam 
axis in a frame in which the pseudorapidities of the two leptons obey 77!+ = —77?-, and 
N(Ei- > Eg+) denotes the number of events for which l~ has a larger energy than £ + . 

This choice of observables is guided by their role in determining different particle proper- 
ties. The first observable in eq. ()2]) is useful for mass measurement, eqs. (J3j) — ([5]) are sensitive 
to the spins of the new particles, and eq. provides information about their couplings. 



3.1 Mass determination 

The variable Mp2 has been proposed for the measurement of particle masses in events with 
two or more invisible objects in the final state [TJ. Mr2 and similar variables have been 
studied extensively in the literature [2], and it was shown that in favorable circumstances 
one can use these variables to determine both the parent mass my as well as the mass of the 
invisible child mx, in particular by including information about initial-state radiation [3]. In 
this paper, therefore, mass determination will not be discussed any further, and the reader 
is referred to Refs. [IHS] for more details. 



3.2 Spin determination 

A useful observable for determining the spin sy of the Y particles is cos 6^ = tanh(A?7«/2), 
see eq. ()3]), which was introduced by Barr in Ref. [5]. It is based on the observation that the 
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final state leptons £ tend to go in the same direction as their parent particles Y^, since on 
average the Y ± are produced with a sizable boost if m Y <C \/s. As a result, the distribution 
of the lepton polar angle 9* Ul in the frame where the pseudorapidities of the two leptons are 
equal in magnitude, is strongly correlated to the production angle 6* between one of the Y 
and the beam axis in the center-of-mass frame. 

The 9* distribution is closely connected to sy. For the spin-0 and spin-1/2 cases one 
finds a characteristic difference which is immediately visible in the formulas 

scalar Y (spin 0): — oc 1 — cos 2 6**, (7) 

fermion Y (spin \): - — n oc 2 + /3 2 (cos 2 6* - 1), (8) 

where (3y is the velocity of the produced Y particles. For spin-1 pair production the situation 
is more complex since here one necessarily needs to take into account a new particle Q in 
the t-channel. Depending on its mass trig, the observable 9* u distribution can be similar 
to the spin-0 case or to the spin-1/2 case, or different from both, as can be seen from the 
numerical results shown in section 13.41 Therefore, in general, the observable ([3]) alone does 
not unambiguously distinguish spin-1 from spin-0 or spin-1/2. 

One advantage of the definition ([3]) is that it is invariant under longitudinal boost, i. e. 
the value of tanh(A^/2) does not depend on the momentum fractions carried by the quark 
and anti-quark in the collision. 

Here we propose two other observables for the determination of the Y spin: the effec- 
tive mass M e ff and the difference between the azimuthal angles of the leptons, A0«, see 
eqs. (jl]) and ([!]). The connection between these variables and sy can be understood from 
the threshold behavior of the partonic cross section qq — > Y + Y~ . If Y + and Y~ are scalars 
they are produced in a p-wave and the cross section behaves like a ~ (3 Y near threshold. 
For fermionic Y, instead, the cross section grows faster near threshold, a ~ (3y. Therefore 
the cross section for fermionic Y pair production reaches its maximum at lower values of 
the Y + Y~ invariant mass, myy, than the cross section for scalar Y pair production. The 
effective mass M e g is strongly correlated to the V-pair invariant mass, and thus the M e g 
distribution will peak at larger values for fermionic Y than for scalar Y (assuming that my 
is equal in both cases and known from measuring the Mt2 distribution). 

The dependence of the cross section on myy also leaves a characteristic imprint on the 
A<pu distribution. Scalar Y pairs will on average be produced with a larger boost than 
fermionic Y pairs. This leads to a more pronounced peak at A0« ~ ir in the scalar case, 
since the larger boost is more likely to produce a back-to-back configuration for the final-state 
leptons, see Fig. [2] 

If Y ± are vector particles, the M e g and A0« distributions depend on the mass ttlq of the 
particle in the t-channel. For itlq ^> my, the Y + Y~ pair production cross section reaches its 
maximum at larger values of myy than both the spin-0 and spin-1/2 cases, since the s-channel 
contribution alone grows monotonically with the center-of-mass energy. In this case, the M e ^ 
distribution for vector Y particles will peak at larger values than the other two cases, and the 
A<pu distribution will be very strongly peaked at tt. On the other hand, Fig. [2] shows that 
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Figure 2: Distributions for M e g (left) and Acfta (right), for combinations 3, 10, and 11 in 
Tab. [H which correspond to Y particles with spin 0, \, and 1, respectively. The plot is based 
on 35000 parton-level events for each combinations without cuts and detector effects, and 
using the mass values my = 300 GeV and mx =100 GeV. For combination 11, two choices 



of the mass of the t-channel particles are shown, 



m 



Q 



1000 GeV and 500 GeV. 



for vtlq of the same order as my the M e g distribution can be similar to either the spin-0 or 
spin- 1/2 cases, depending on the precise value of rriQ. Nevertheless, even for a relatively low 
value rug = 500 GeV@ the A(j)u distribution is still distinctly different for spin-1 compared 
to the other spin cases. By using all three observables (J3I) — (J5]) in combination one therefore 
obtains the best discrimination power and can unambiguously distinguish between sy = 0, 
2, and 1. 

Similar to tanh(Ar]«/2), also A0« and M e s are invariant under longitudinal boosts, and 
thus very well suited for hadron colliders. 

It needs to be pointed out that the three variables, tanh(A?7^/2), A0« and M e g, are 
primarily sensitive to the spin of the parent particle Y, but not of the child particle X. 
Indeed, as can be seen from the numerical results in sections 13.41 and HJ it is very difficult to 
independently determine the X spin. 

3.3 Coupling determination 

Experiments at LEP and SLC have determined the couplings of the Z boson to SM fermions 
with very high precision, in particular by measuring various left-right and forward-backward 
asymmetries |15j . 

Similarly, for the class of processes corresponding to Fig. [TJ one can in principle try to 
extract information about the ZYY coupling by constructing a forward-backward asymmetry 
for pp — > Y + Y~ at the LHC Although the initial pp state is symmetric, the incoming quark 
for a gg-initiated process often stems from one of the valence quarks of the protons and 
thus tends to have a larger momentum than the incoming anti-quark. Therefore one can 
define the forward direction by the direction of the overall longitudinal boost of an event. 

$For such low values of toa one should see a signal from direct production of the Q particle at the LHC. 
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Form of the XY£ coupling: 
• YX^i 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

^ZA 

Figure 3: Relationship between the coupling ratio Rza = g{ZYY)/e and the asymmetry 
Af+p-, for two different chiralities of the interaction in the Y ± — > £ ± X° decay, where Y is a 
fermion and X is a vector boson. The combinations 6 and 9 in Tab. [JJ are indicated by the 
open square and circle, respectively. The plot is based on parton-level results without cuts 
and detector effects, and using the mass values my = 300 GeV and mx = 100 GeV. 

However, since we neglect effects from electroweak symmetry breaking in the new physics 
sector, all combinations in Tab.UJhave parity-even ZYY couplings and the forward-backward 
asymmetry for pp — > Y + Y~ is exactly zero. 

On the other hand, the coupling between the incident qq pair and Z boson has a parity- 
odd axial-vector part, which results in the Y + Y~ pair being produced with a non-vanishing 
polarization asymmetry (unless Y ± are scalars). This polarization asymmetry can be probed 
through the decay Y ± — > i^X , since the interaction responsible for the decay is either left- 
or right-handed and thus sensitive to the Y polarization, see Tab. [TJ 

In the center-of-mass frame of the Y + Y~ system this leads to a forward-backward asym- 
metry for the final-state leptons. As mentioned above, in the lab frame the forward direction 
is defined by the overall boost of an event, which is closely correlated to the direction of the 
more energetic of the two leptons. Therefore we define the observable given in eq. ([6]), 

_ N(E e - > Eg.+ ) — N(E e + > E e -) 

~ N(E e - > E e +) + N(E e + > E e -)' m 

The asymmetry is partially washed out by the mass my, which can cause a spin flip before 
the Y decays, but we expect a non-vanishing result as long as m Y y/s. 

Eq. [6] is mostly useful for discriminating between combinations with Sy = |, since scalars 
do not carry any polarization and lead to a vanishing asymmetry, and there is only one 
combination with vector Y particles in Tab. [TJ The value of A^+i- is connected to the size of 
the ratio Rza = g(ZYY)/e between the ZYY and ^YY couplings and to the sign of the 75 
term in the XY£ coupling. This is illustrated in Fig. [3] for the case that Y is a fermion and 
X is a vector boson. As evident from the plot, there is a strong correlation between Rza 
and Ai+£-, but one can encounter a two- to three-fold ambiguity when trying to determine 
Rza from the measured value of A1+1-. 
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Table 3: yx^ values for a 5-bin x 2 -test to discriminate between pairs of model combinations 
with different spin of the parent Y particle. The combinations 3, 10, and 11 from Tab. [JJ have 
been chosen as examples of models with Y particles of spin 0, |, and 1, respectively. Model B 
is assumed to represent the simulated "data", while model A is the test hypothesis. The 
results are based on samples of 5000 parton-level events without cuts and detector effects, 
and using the input values my = 300 GeV, mx = 100 GeV, and \fs = 14 TeV. 

3.4 Numerical results 

Using CompHEP we have generated parton-level events for all 11 combinations in Tab. [TJ 
for a center-of-mass energy of yfs = 14 TeV and my = 300 GeV and mx = 100 GeV. In the 
following discussion we will assume that the Y and X masses are known from observables 
like Mt2, and for simplicity the uncertainty in the mass determination will be neglected. 

As explained in the previous section, the observables tanh(A?7«/2), M e g, and A(j)gg, see 
eqs. (JHJ) (J5j) can be used to determine the spin sy of the parent particle Y. We have checked 
this by performing a 5-bin x 2 analysis for a sample of 5000 events for each model combi- 
nation, assuming Poisson statistics for the statistical error. Table [3] shows the results for 
combinations 3 (with scalar Y), 10 (with fermion Y), and 11 (with vector Y) as examples. 
For the case of vector Y particles, results for two sample values of the t-channel fermion 
mass mg are given. As evident from the table, by combining the three observables, one can 
distinguish all three spin combinations from each other with a significance of more than 18 
standard deviations, for the given choice of masses and total event count. This is true even 
for relatively small values of m^ ~ (9(0.5 TeV). 

We have checked that the results are very similar if combinations 3 or 10 are replaced by 
any of the other combinations with spin-0 or spin-1/2 Y particles, respectively. Furthermore, 
for any two models with identical sy the distributions for all three variables are statistically 
consistent, irrespective of the spin of X. 

To discriminate between models with identical sy but different SU(2) representations 
of the Y and X particles one can take advantage of the charge asymmetry Ag+g- in (Q. 
As mentioned in the previous section, one cannot obtain a non-zero asymmetry if the Y 
particles are scalars, and we have confirmed this statement explicitly with our simulation 
results. However, for sy = \, Ag+g- can yield useful information about the structure of 
the ZYY and XYl couplings. Results for the total asymmetry, without cuts or acceptance 
effects, are given in Tab. H]for all combinations with fermionic Y in Tab. [TJ 
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Table 4: Values for the asymmetry At+%- for combinations with fermionic Y in Tab. [I] based 
on simulated parton-level events for m Y = 300 GeV, m x = 100 GeV, and y/s = 14 TeV. 
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Table 5: Statistical signihcance, in units of standard deviations, for the discrimina- 
tion between combinations with fermionic Y in Tab. [I] using the differential asymmetry 
dAi+g- jd tanh( Arju/ '2) . Numbers in bold face indicate a difference of at least 20 standard 
deviations, while gray italic numbers denote a significance of less than three standard devia- 
tions. The results are based on samples of 5000 parton-level events without cuts and detector 
effects, and using the input values m Y = 300 GeV, m x = 100 GeV, and y/s = 14 TeV. 

In general, Ai+i- becomes maximal for events with large values of cos 6^ = tanh(Ar/ a /2) 
close to 1, i. e. when the Y + Y~ pair is produced in the forward/backward direction. However, 
this correlation between A t + t - and tanh(Ar/^/2) depends to a lesser extent also on spin 
effects in the decay Y — > IX and thus can be markedly different for models with opposite 
chirality of the XYl vertex. As a result, the significance for distinction between such models 
is increased by performing a binned analysis for the distribution dAi+i-/dta,nh(Ar}u/2). It 
turns out that the highest sensitivity is obtained by using just two bins. 

Table lists the statistical significances for discriminating between any pair of the com- 
binations 4-10 from Tab. [1] based on this observable. Models that have different signs for 
Ai+£- can be distinguished with more than 20 standard deviations for an signal event sample 
of 5000 events (bold face numbers in the table). 

However, the combinations 4, 7, and 10, as well as 6 and 9, are indistinguishable at the 
three-sigma level (gray italic numbers in the table). It turns out that also when considering 
any other variables in eqs. ([2])—© one cannot achieve a higher significance for discriminating 
between these models. 

Note that the variable Ai+i- has some sensitivity to distinguish between models which 
differ only through the spin of the X particle, i. e. between combinations 4 and 6, or 7 and 9 
in Tab. [TJ Assuming a signal sample of 5000 events, as in Tab. \5\ a discrimination significance 
of about five standard deviations can be achieved for these pairs. 
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(model A, mo 


del B) 




VcLIldDlt3 


{6,l\J) 


(3,H) 
[M d =l TeV] 


(3,11) 
[M a =0.5 TeV] 


(10,11) 
[M 6 =l TeV] 


(10,11) 
[M a =0.5 TeV] 


tanh(Ar/«/2) 


23.2 


20.1 


30.6 


6.3 


10.1 


M eff 


40.0 


9.7 


24.2 


36.8 


12.6 




27.2 


15.5 


6.0 


40.7 


23.8 


All combined 


40.0 


20.1 


30.6 


40.7 


23.8 



Table 6: \J~x* values for a 5-bin x 2 -test to discriminate between pairs of model combinations 
with different spin of the parent Y particle, for a sample of 5000 events passing the detector 
simulation and selection cuts in eq. fl9]). The notation and input parameters are the same as 
in Tab. 



3.5 Simulation results 

The analysis in the previous section does not take into account detector acceptance and 
signal selection cuts. To study the influence of these effects on the results we have passed the 
parton-level events generated by CompHEP [TJ] through Pythia 6.4 [16] and PGS4 [T7] . 
By including initial-state radiation and parton showering in the Pythia simulation one can 
furthermore evaluate whether fluctuations of the initial-state transverse momentum might 
wash out the characteristic features for the model discrimination. 
In Ref. [5] it has been shown that the selection cuts 

N(£+) = N(r) = 1, met > 150 GeV, max{p T/ ±} > 40 GeV, min{p T/ ±} > 30 GeV, 
^ T > 100 GeV, M T2 > 100 GeV, |pV + p T/+ + p T/ -| < 100 GeV, 
Ptj < 100 GeV, N b = 0, (9) 

reduce the SM background rate to about 1.6 fb. Here A^(£ ± ) denotes the number of visible 
leptons £ ± = e ± , /x ± in the central detector, Nb denotes the number of vertex b tags, ma is the 
di-lepton invariant mass, and px,j refers to the transverse momentum of any reconstructed 
jet. With these cuts one obtains a selection efficiency for the signal process pp — > Y + Y~ — > 
£ + £~XX between 27% and 40%, depending on the specific type of Y and X particle. As 
listed in the appendix, this corresponds to measurable signal cross sections between about 
1 fb and 200 fb. For concreteness we will assume 5000 observed events, which corresponds 
to the expected yield of model combination 7 for an integrated luminosity of 200 fb -1 . In 
comparison, the SM background of about 300 events is small and can be neglected. 

Tables M arid [7] summarize the significance for distinguishing between models with differ- 
ent Y spin and with different couplings, assuming 5000 measured events for y/s = 14 TeV, 
my = 300 GeV and mx = 100 GeV. Overall, the obtained significances for the spin dis- 
crimination are comparable to the parton-level results in Tab. [31 and in a few cases the 
significance is even higher. This seemingly surprising outcome is related to the fact that we 
compare the same number of "observed" events in the previous section and in this section, 
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9 




5 


23 
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3.3 


20 










PQ 


7 


2.0 


22 


1.6 
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8 


22 


2.1 


19 


21 








9 


3.7 


19 


1.7 


3.3 


17 






10 


1.3 


25 


4.1 


2.1 


23 


5.3 



Table 7: Statistical significance, in units of standard deviations, for the discrimination 
between combinations with fermionic Y using dAg+£- /dtanh(Ar/^/2). The results are based 
on a sample of 5000 events passing the detector simulation and selection cuts in eq. 0, with 
notation and input parameters are the same as in Tab. 

but in the latter case the cuts remove part of the phase space, leaving a higher event yield 
in the remaining phase-space region. 

For the coupling determination one finds that the asymmetry Ai+i- is washed out no- 
ticeably by the cuts, leading to substantially reduced significances in Tab. [7] compared to 
Tab. [51 Nevertheless, models with different sign for Ai+i- can still be distinguished with at 
least 17 standard deviations. 

In summary, for most cases, selection cuts and smearing effects only moderately affect 
the capability for identifying particle properties with the described observables. Of course 
the selection cuts reduce the overall event number, which however also depends on the 
model-dependent total cross section and thus is left as a free parameter here. 

4 Comparison with Automated Likelihood Analysis 

An alternative approach for the analysis of a new-physics signal is an automated likelihood 
test for a sample of measured events. With such a computerized procedure it is in general 
not possible to clearly separate properties like spin and couplings, but it offers the advantage 
of reaching a higher sensitivity by using the complete event information instead of specific 
observables. A very appealing realization of an automated likelihood analysis is the Matrix 
Element Method (MEM) [HUE], which uses parton-level matrix elements to specify the 
theoretical model that is compared with the data. The method can be used to measure one 
or several parameters of the model by finding the maximum of the likelihood for a sample of 
events as a function of these parameters. As of today, the MEM achieves the most precise 
determination of the top-quark mass [TU] and new-physics particle masses [I]. 

For each single event, with observed momenta p™ 3 , the MEM defines a likelihood measure 
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(model A, model B) 


(3,10) 


(3,11) 
[M =l TeV] 


(3,11) 
[M a =0.5 TeV] 


(10,11) 
[M d =l TeV] 


(10,11) 
[M 6 =0.5 TeV] 


60 


59 


61 


85 


87 



Table 8: Statistical significance, in units of standard deviations, for the discrimination 
between pairs of model combinations with different spin of the parent Y particle, based on 
the MEM. A sample of 5000 parton-level events without cuts and detector effects has been 
used. The notation and input parameters are the same as in Tab. 



that it agrees with a model for a given set of model parameters a: 

M a ( Pt )\ 2 l[5(p t -pr). (10) 

iGvis 

Here fx and j% are the parton distribution functions, M a is the theoretical matrix element, 
and a a is the total cross section, computed with the same matrix element. The 3- momenta 
pj" s of the visible measured objects are matched with the corresponding momenta p { of the 
final state particles in the matrix element, while the momenta of invisible particles (weakly 
interacting particles, such as the X particle in our case) are integrated over. 

For a sample of iV events, the combined likelihood is usually stated in terms of its 
logarithm, which in the large- N limit can be interpreted as a x 2 value, 

N 

X 2 = -21n(£) = -2 5>P(P?il«), (11) 

n=l 

where are the measured momenta of the nth event. 

The MEM is particularly useful for signals that cannot be fully reconstructed due to 
invisible final-state particles, and it can be applied to determine the masses of both X and 
Y in processes of the type in eq. ([Tj) [3j. Here we will assume that the masses are already 
known and instead focus on the discrimination between the models in Tab. [TJ 

Matrix elements for all 11 combinations in the table have been computed with the help of 
CompHEP and implemented into a private code for performing the phase-space integration 
in ( II Op . Similar to section 13.41 only parton-level events without cuts have been used in this 
analysis. Results for model comparisons are listed in Tables M and |9j 

As can be seen from Tab. EJ the MEM achieves a much higher significance for discriminat- 
ing between combinations with different sy, see Tab.[3]for comparison. This is not surprising 
since several observables, eqs. ([I])-©, were found to be sensitive to the Y spin, indicating 
that none of them captures all relevant information. Note also that the results in Tab. [8] do 
not depend strongly on the unknown mass of the t-channel fermion Q for combination 11. 

The MEM can also distinguish between combinations that all have spin-1/2 Y parti- 
cles but which differ in the SU(2) representations of X and Y, as shown in Tab. IHJ It is 
interesting to note that in most cases the statistical significance achieved by the MEM is 
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8.3 


25 
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2.5 


29 


9.1 
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8 


30 


2.3 


27 


28 






a 


9 


8.6 


26 


3.0 


9.0 


27 






10 


15 


41 


22 


14 


43 


20 



Table 9: Statistical significance, in units of standard deviations, for the discrimination 
between combinations with fermionic Y in Tab. [I] based on the MEM. A sample of 5000 
parton-level events without cuts and detector effects has been used. The notation and input 
parameters are the same as in Tab. 

comparable to the results obtained with the asymmetry Ag+i- in Tab. [5j An exception is 
combination 10 which can be distinguished from the other combinations with substantially 
higher significance using the MEM compared to Ap+i-. This implies that the asymmetry 
Ai+ k - captures essentially all measurable information about the ZYY and XY£ couplings, 
except for the special case of model 10. 

Similar to the results of the previous section, it is found that one cannot discriminate 
very well between combinations with Y singlets and Y doublets, i. e. between 4 and 7, 5 and 
8, or 6 and Likewise, the MEM results for the combinations 1, 2, and 3 with scalar Y 
differ by less than one standard deviation, and thus are completely indistinguishable. 

5 Conclusions 

This paper presents a comprehensive analysis of new physics processes of the form pp — > 
Y + Y~ — > £ + £~X°X°(£ = e,fi), where X° is stable and weakly interacting, leading to a 
signature of two opposite-sign same-flavor leptons and missing momentum. To minimize 
model assumptions, all possible combinations for the spins and weak SU(2) couplings of X 
and Y have been considered, allowing for spin 0, | and 1, and SU(2) iso-singlets, -doublets 
and -triplets, see Tab. [TJ 

The signal processes have been analyzed with two different and complementary ap- 
proaches. The first method is based on specific observables. Concretely, we have studied 
three variables for the measurement of the spins and one asymmetry for the extraction of in- 
formation about the couplings of the new particles. Secondly, an automated strategy called 
the Matrix Element Method has been used, which algorithmically computes a likelihood 

''Note, however, that a better differentiation between these cases would in principle be possible with more 
statistics, requiring significantly larger amounts of integrated luminosity. 
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that a given event sample agrees with some model interpretation supplied in the form of a 
theoretically calculated matrix element. 

It has been found that the spin sy of the parent particle Y can be determined with 
high statistical significance, so that a sample of a few hundred signal events is sufficient for 
discrimination at the 5cx-level. Furthermore, it was shown that the asymmetry A^+p- defined 
in eq. is instrumental in distinguishing between model combinations that all have sy — \ 
but different Y and X couplings. The majority of possible coupling assignments can be 
differentiated with high significance, but it turns out that for the pairs 4 and 7, 5 and 8, as 
well as 6 and 9 in Tab. [I] one cannot achieve a 3a discrimination with a realistic number of 
a few thousand events. This is related to the fact that the relationship between the ZYY 
coupling strength and the observable asymmetry is not monotonic and can involve degenerate 
solutions. Remarkably, the same model combinations are also difficult to distinguish with the 
Matrix Element Method, which demonstrates that the asymmetry Ap+p~ reflects all relevant 
information about the couplings of the underlying model. 

For sy = it is generally impossible to discriminate between cases with different couplings 
or with different spin of the X particles, due to the absence of spin correlations between the 
production and decay stages of the process. For sy = 1 the coupling structure of the process 
is essentially fixed by gauge invariance and thus already uniquely known once the vector 
nature of Y has been determined. 

Our findings indicate that even for the challenging case of a process with a short one-step 
decay chain it is in general possible to separately determine the spins and couplings of the 
new heavy particles. The results in this paper have been presented for the specific choice 
of masses my = 300 GeV and mx = 100 GeV, but we have checked explicitly that the 
essential features are unchanged for my = 200 GeV. While the main goal of this study was 
the development of the theoretical framework and conceptual ideas, we have also performed 
a fast detector simulation with selection cuts for the suppression of standard model back- 
grounds and found that qualitatively our conclusions still hold. Nevertheless, a dedicated 
experimental simulation with a careful evaluation of systematic errors, including the influ- 
ence of uncertainties in the Y and X masses, would be required to check the viability of our 
results under realistic conditions. 
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Appendix: Model cross sections 

The following table lists the tree-level parton-level production cross sections a pio d for the 
process pp — > Y + Y~ — > £ + £~X°X°(£ = e,fi), for the 11 independent combinations from 
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Tab. [TJ Also shown are the measurable cross sections <7 meas after inclusion of detector effects 
and the cuts in eq. The cross sections have been computed with CompHEP. 



Combination 


Oprod [fb] 




1 


3.62 


1.45 


2 


8.50 


3.36 


3 


9.65 


3.11 


4 


41.4 


11.45 


5 


41.4 


11.70 


6 


41.4 


14.05 


7 


89.6 


25.0 


8 


29.9 


8.47 


9 


89.6 


31.4 


10 


112 


31.2 


11 [Mq=0.5 TeV] 


179 


48.3 


11 [M =l TeV] 


445 


137 
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